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We simulate the kinetic Ashkin- Teller model with both ordered and 
^ ■ disordered initial states, evolving in contact with a heat-bath at the 

critical temperature. The power law scaling behaviour for the magnetic 
order and electric order are observed in the early time stage. The 
values of the critical exponent 6 vary along the critical line. Another 
dynamical exponent z is also obtained in the process. 
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It has been predicted by analytical calculations |^ and supported in Monte 
Carlo simulations 0, ^ |, H that there exists scaling in the macroscopic short- 
time regime for some critical dynamic processes. The universality was confirmed 
for the Ising model and the Potts model with various dynamics and different 
lattice structures 0, H- Basing upon the universal scaling hypothesis for the 
initial stage, promising methods have been proposed for numerical measurements 
of critical exponents, including both static and dynamic ones [P, |10|, |TlI]. It has 



been suggested that the critical point can be also determined in this stage [11]. 
The universal behaviour of the short-time dynamics is found to be quite general 



12, 13, 14, 15], e.g., in connection with ordering dynamics or damage spreading 



1^, [T^ and surface critical phenomena [|l| 



It is interesting to see the short-time behaviour of models possessing contin- 
uously varying critical exponents in equilibrium. It was found by Kadanoff and 
Wegner that there is a connection between the continuous variation of the critical 



exponents and the existence of a marginal operator ]jT9[. The variation of critical 
exponents in the early time evolution implies that the operator keeps in marginal 
even in the non-equilibrium initial stage. 

The Ashkin- Teller model is one of such models that has week universality, i.e., 
the critical exponents vary with parameters of the interaction ]^. It contains 



two species of spins {cTj = ±1} and {rj = ±1}, located on a square lattice. The 
interaction is described by the Hamiltonian 

H = J2 l^i^i^j + ^i^j) + K^aiajTiTj)] (1) 

<i,j> 

This model is dual to the solved Baxter's eight vertex model. It has five phases. 
We only investigate the exactly known segment of the critical line that separates 
the disordered phase and the phase where both the magnetic order and electric 
order are not zero. Along that line, the critical exponents vary continuously. The 
critical line is given by the equation 

exp(-2K4) = sinh(2if) (2) 

where (ln3)/4 < i^' < -|-cxo. Besides the critical exponent u that governs the crit- 
ical behaviour of the correlation length, there are other two independent critical 
exponents jSm and /5e corresponding to two order parameters, the magnetic order 
< (jj >, and the electric order < CiTi >, respectively. The critical exponents for 
the equilibrium state have been known as [^, |2^ 



/3™ = ^. A = T^ (3) 



3-2?/' 24-16?/' 12-8?/ 

where the parameter y, ranging from to 4/3, is related to K4 by the equation 

cos(7r?//2) = i[exp(4ir4)-l] (4) 
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When dynamics is introduced into the Ashkin- Teller model, besides the well- 
known z that describes the divergence behaviour of the time-correlation length, 
there are two critical exponents related to the dimensions of two initial orders. 
As argued by Janssen et al., the initial scaling emerges at fixed point under 
the renormalizational group transformation. By naive coarsing transformation, 
one can see two fixed points for initial states, i.e., mo = Cq = at very high 
temperature and mo = Cq = 1. Our simulations are around these two fixed 
points. 

The magnetic moment M^^^ and the electric moment E^^^ are defined as 

=< [(5: a,)lL^f >t, E^'^ =< [(^ a,T,)/L']' >, (5) 

i i 

In analogy to the Ising model, we suppose the following finite-size scaling relations 
hold in the vicinity of the first fixed point, 

M^^\t,L,mo) =h-^^^/''M^^\h-'^t,h-^L,h''^mQ) (6) 

for Co = 0, and mg small, and 

E^^\t, L, eo) = b~'"^'^/''E^''\b-'H, b'^L, fo^-^co) (7) 

for mo = 0, and eo small. We have denoted the anomal dimensions of the initial 
magnetization and electric order as Xm and Xe respectively. For k = 1, taking 
b = t^/' and assuming the initial orders are small enough, we can expand the 
magnetization(electric order) with respect to mo(eo) and obtain 

M{t) = mQt^-^Fra{t-^'"^L) (8) 

E{t) = eot'^F,it-'/'^L) (9) 

where 9m = — Pm/ ^) / Zm, and 9^ = (xe — Pe/ ^) I -2e- The scaling functions 
and Ff. take account of the deviation from power laws due to the finite-size effect. 
They tend to constants as L tends to infinity and time tends to zero, but bigger 
than the microscopic time scale tmic- It has been shown in various models that 
traic is ignorable in the heat-bath dynamic process. In the following fits, we will 
consider F^ and F^, as constants. 

It has been stressed in Ref. that the initial states must have very short 
correlation lengths and that a sharp preparation of the initial state improves the 
result. For the Ashkin- Teller model, there are two initial order parameters, mo 
and Co. So one has more than one way to approach the disordered fixed point. For 
instance, we can either let eo = and mo small or reverse, mo = and eo small, 
later on referred to as initial condition I and initial condition II, respectively. 
We find that the relaxation patterns from these two initial conditions are very 
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L = 180 


L = 60 


L = 30 


y 


mo = 0.02 


mo = 0.02 


mo = 0.04 


mo = 0.06 


mo = 0.08 


mo = 0.04 





-.020(1) 


-.021(4) 


-.022(5) 


-.019(5) 


-.021(2) 


-.022(6) 


1/6 


-.013(1) 


-.014(2) 


-.014(5) 


-.013(4) 


-.014(2) 


-.015(6) 


1/3 


.009(2) 


.008(1) 


.007(5) 


.007(3) 


.007(3) 


.007(4) 


1/2 


.041(5) 


.041(3) 


.041(5) 


.041(2) 


.041(2) 


.039(6) 


2/3 


.085(4) 


.085(4) 


.085(5) 


.084(3) 


.083(2) 


.086(6) 


5/6 


.137(2) 


.137(2) 


.134(3) 


.136(3) 


.134(1) 


.135(3) 


1 


.191(2) 


.189(2) 


.188(3) 


.188(2) 


.185(1) 


.187(3) 


7/6 


.230(1) 


.229(2) 


.227(1) 


.225(1) 


.223(1) 


.225(3) 



Table 1: Results for 6m with initial condition I of mo = 0.02 on L = 180, and of 
mo = 0.02, 0.04, 0.06, and 0.08 on L = 60, mo = 0.04 on L = 30. In each case, 
the magnetization in the time interval 1 < t < 5 is used for the least-square fit 
to the power law. 

different. However, for both initial conditions, power laws are observed. The 
exponents 6 depend on the initial conditions and vary with y. 

The short-time behaviour of the magnetization with the initial condition I 
can be seen in Fig. 1. The points are magnetizations averaged over 40000 inde- 
pendent samples with L = 180 and mo = 0.02. Since the two species of spins 
are symmetry, the effective sample number is doubled. Samples of each point are 
grouped into 4 runs, and the errors are obtained from them. The parameter y 
is increasing from the bottom to the top. The lines are curves of the power law 
Eq. (P) with the best-fit exponents 6m to the magnetizations in the first five time 
steps as given in Tab. 1. We stop at t = 5 for that gives the smallest fluctuation. 
Tab. 1 also contains 6^, for L = 30, and 60, which shows that the finite-size effects 
are smaller than the statistical fluctuation. By comparing results for various mo 
from 0.02 to 0.08, on lattice L = 60, one can see that the 6m are quite stable. At 
the decoupled point y = 1, the Ashkin- Teller model reduces into two independent 
Ising models. Our best value 6m = 0.191(2) should be compared with the existing 
numerical results for the Ising model from Refs. 0, [l^, 6 = 0.191(1), 0.191(3), 
and from those obtained from auto-correlation before ^. 

Figure 2 shows the initial relaxation of the electric order starting from the 
initial condition II with eo = 0.02, L = 180 in the double-log scale. Each point 
is an average over not less than 80000 independent samples. To obtain stable 
results for y = 1, we used up to 480000 independent samples. The lines from 
the top to the bottom are curves of the power law Eq.(^ with the best-flt 6e, 
as given in Tab. 2, corresponding to y = 0,1/6, 1, respectively. For y < I, 
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L = 180 


L = 60 


L = 30 


y 


Co = 0.02 


Co = 0.02 


Co = 0.04 





-.021(6) 


-.022(5) 


-.021(3) 


1/6 


-.036(6) 


-.037(6) 


-.036(3) 


1/3 


-.082(7) 


-.085(5) 


-.083(3) 


1/2 


-.169(6) 


-.172(6) 


-.171(3) 


2/3 


-.308(9) 


-.312(11) 


-.306(4) 


5/6 


-.510(9) 


-.504(10) 


-.498(5) 


1 


-.812(30) 


-.793(29) 


-.760(13) 



Table 2: Results for 6^ with initial condition II of cq = 0.02 on L = 60 and 180, 
and of Co = 0.04 on L = 30. The electric order in the time interval 1 < t < 5 is 
used for the least-square fit to the power law. 

the order monotonously decreases with time in a power law. Discontinuity is 
found around y = 1, which is the turning point of coupling K^^ from positive 
to negative. For y = 7/6, the order jumps to a negative value at the first time 
step. The fiuctuation for the electric order is much bigger than that for the 
magnetization since the electric order decays rapidly to very small values. That 
is the reason why we have to measure the exponent 9e at the very beginning of 
the relaxation (1 < t < 5). 

Now, we turn to the measurement of the dynamic exponent z. Traditionally, 
z is measured from the long-time exponential decay of the time correlation or the 
magnetization of the systems [0, Due to the critical slowing down, this is 
somehow difficult. Recently, a few methods have been proposed to estimate z in 
the short time stage of dynamic processes. Stauffer suggested to obtain z from 
the power law decay of the magnetization with the known static exponent (3/v 
as input. For a large enough lattice, one may expect a power law decay of the 
orders 

M{t) ~ t-^™/'^^- (10) 
E{t) ~ t"/^-/'^^- (11) 

before the exponential decay starts. For the Ising model and the Potts model with 
ordered initial states, it has been confirmed that the power law decay happens in 
the short-time regime [^, |2^, |10|, |TT1. The exponent z may also be measured in- 
dependently from the time- dependent Binder cumulants in the short-time regime 
with either ordered or disordered initial states as proposed in Refs. |P, |lOl . Since 



higher order moments are needed, the results are relatively more fiuctuating. 
With the initial dynamic exponent 9 in hand, Ref. suggested to determine z 
from the power law behaviour of the autocorrelation. In this paper we do not 
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use this method since the 9e for the Ashkin- Teller model is big and difficult to be 
precisely evaluated. 





L = 64 


L = 90 


y 


Zm 


Ze 


Zm 


Ze 





2.23(4) 


2.24(4) 


2.21(2) 


2.21(3) 


1/6 


2.27(6) 


2.37(6) 


2.24(5) 


2.34(5) 


1/3 


2.31(2) 


2.37(2) 


2.29(4) 


2.35(5) 


1/2 


2.31(4) 


2.32(4) 


2.26(3) 


2.27(3) 


2/3 


2.26(2) 


2.24(2) 


2.24(4) 


2.23(3) 


5/6 


2.19(2) 


2.17(2) 


2.20(1) 


2.19(1) 


1 


2.17(1) 


2.17(1) 


2.15(2) 


2.15(2) 


7/6 


2.12(3) 


2.14(3) 


2.14(3) 


2.16(4) 



Table 3: Results for z obtained by fitting the power-law decay of the order 
parameters from an ordered initial state with L = 64 and L = 90, where Zm and 
Zf, are corresponding to the magnetization and the electric order respectively. The 
exact values of 13^1 t-' and (3e/i' are used as input. To avoid the possible deviation 
from the power-laws in the beginning of the relaxation, the ffist 100 time steps 
are skipped in the fitting. 

In Figs. 3 and 4, we plot the time evolution with the ordered initial state 
for the magnetization and the electric order respectively, in double-logarithmic 
scale. Each point is measured from 16000 samples on a lattice with L = 90. 
The curves are the power laws as given by Eqs. (|l^) and (pUD with the best fit 
powers for various y. To avoid effects of tmic, the first 100 time steps have been 
skipped. With the static exponents Pm/^ and Pe/j^ given by Eq.(||), we obtain 
the dynamic exponent Zm and Ze from the powers of Fig. 3 and 4 respectively. 
The results are compared with those of L = 64 in Tab. 3. The finite-size effects 
are about 2.5%, comparable with the statistical errors. As y varies from zero to 
7/6, a remarkable change of z may not be accounted only to finite-size effects and 
statistical errors. This may indicate that the exponent z is varying. However, 
within the statistical errors, we can not distinguish between Ze and Zm- At the 
decoupling point y = 1, we obtain the dynamic exponent z = 2.15(2) from both 
magnetic and electric powers, which is consistent with recent result for the Ising 
model obtained by other authors [^, |2^, |2^, |28|, 

We also estimate the dynamic exponent z independent of the static exponents 
by the finite-size scaling fit of the time- dependent cumulants on lattices with 
L = 64 and L = 90. The method has been described in detail in Ref. |rD[. We 



get the mean value z = 2.16(6) from the global fit and z = 2.16(8) from the local 
fit. The large fiuctuation prevents us to make any conclusion on whether z is 
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varying. Using the mean value z — 2.16 as input, the static exponents Pm/j^ and 
Pe/j^ can be estimated from the finite-size scaling of the magnetization and the 
electric order respectively. The results turn out to be consistent with the analytic 
results with deviation smaller than 10%. 

In conclusion, we have confirmed the initial scaling for the kinetic Ashkin- 
Teller model which is known to possess weak universality in the equilibrium. The 
power laws fit to the measurements remarkably well in the short-time evolution 
of the order parameters with both disordered and ordered initial states. The 
exponents 9m and 9^ are found varying with y. The initial order increase is only 
observed for y > 1/6 with the initial condition 1. This fact can not be explained 
by the straightforward mean-field argument that claims the initial order increase 
should happen when the critical temperature is lower than the mean-field critical 
temperature. It seems that more careful analysis is needed. 

For bigger y{> 1/6), Cq has negative dimension Xg. This means that the time 
scale associated with cq is a short-time scale, in contrast with the long-time scale 
associated with mo. In Fig. 2 one can see that the scaling window gets narrower 
and narrower as y increases (the dimension decreases). When the time scale of 
Co is comparable with the microscopic time scale, one will not see the initial power 
law behaviour of the electric order. On the other hand, for y = the electric 
order and the magnetization have the same value for 9 within the fiuctuations. 
This should be clear since at that point K — K4, i.e., there is symmetry between 
two orders. 

The data in Tab. 3 indicate that the exponent Zm is approximately equal to 
Ze- This may imply that the model has only one time correlation length which 
goes to diverge at the critical temperature, as there is only one space correlation 
length in the equilibrium Ashkin- Teller model as indicated by the equality of 
for two orders. 

The exponent z may vary with y, however, the errors are still too big to allow 
a definite conclusion. There are still some interesting questions left for further 
study. For instances, one should better understandthe infiuence of the choic of 
the initial conditions, the discontinuity of the electric order at the decoupling 
point, the tricritical point y = 0, etc. To locate the critical line segments which 
are not exactly known would be also a good way to test the power of short-time 
dynamics. 
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Figure captions 



1. Magnetization versus time for the initial condition I with mo = 0.02, cq = 0, 
and L — 180. The curves from the bottom to the top are best fits to 
the power law Eq. (8) with y = 0, 1/6, 1/3, 1/2, 2/3, 5/6, 1, and 7/6, 
respectively. 

2. Electric order versus time for the initial condition II with mo = 0, Co = 0.02, 
and L = 180. The curves from the top to the bottom are best fits to the 
power law Eq. (9) with y = 0, 1/6, 1/3, 1/2, 2/3, 5/6, and 1, respectively. 

3. The magnetization for various y on L = 90, left for relaxation from the 
ordered state. The curves are best fits to the power law Eq. (10) with Zm 
given in Tab. 3 and Pm/j^ given by Eq. (3) as input. 

4. The electric order for various y on L = 90, left for relaxation from the 
ordered state. The curves are best fits to the power law Eq. (11) with Ze 
given in Tab. 3 and /3e/i^ given by Eq. (3) as input. 
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Figure 1 
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Figure 2 
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Figure 4 
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